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EQUATIONS  FOR  DETERMINING  THE  GAS 
AND  BRINE  VOLUMES  IN  SEA  ICE  SAMPLES 


Cordon  F.N.  Cox  and  Wilford  F.  Weeks 


INTRODUCTION 

A  knowledge  of  the  total  porosity  of  sea  ice  is 
important  in  interpreting  its  physical  properties.  It 
has  been  shown,  both  theoretically  and  experiment¬ 
ally,  that  the  mechanical,  thermal,  and  electrical 
properties  of  sea  ice  vary  with  the  relative  volume 
of  brine  existing  in  the  ice  (Weeks  and  Assur  1 967, 
Schwarz  and  Weeks  1 977).  However,  the  air  or  gas 
volume  in  the  ice  is  also  important.  This  is  particu¬ 
larly  true  for  low  salinity  sea  ice,  such  as  multiyear 
ice,  where  the  gas  volume  may  make  up  a  major  por¬ 
tion  of  its  total  porosity. 

The  gas  volume  of  sea  ice  can  be  calculated  given 
the  density,  salinity,  and  temperature  of  the  ice. 
Generally,  Assur’s  (1958)  phase  equilibrium  table 
(Table  III,  p.  124)  is  used  to  perform  such  calcula¬ 
tions.  This  table  gives  the  relative  masses  of  the  dif¬ 
ferent  components  in  "standard  sea  ice”  at  tempera¬ 
ture  intervals  of  two  degrees  Celsius.  These  data  are 
primarily  based  on  the  experimental  determinations 
of  Nelson  and  Thompson  (1954)  and  Ringer  (1928). 
Here  the  expression  “standard  sea  ice"  is  used  to  refer 
to  sea  ice  of  such  a  composition  that  the  relative 
concentrations  of  ions  in  its  meltwater  are  the  same 
as  in  normal  sea  water.  From  the  density  of  the 
various  components,  the  theoretical,  or  gas-free  den¬ 
sity  of  the  sea  ice  can  then  be  calculated  and  com¬ 
pared  to  the  measured  density  to  determine  the 
actual  gas  volume  present  in  the  ice.  However,  since 
the  table  is  based  on  standard  sea  ice  having  a  salinity 
of  34.325  °/oo,  it  is  first  necessary  to  reduce  the  data 
to  the  appropriate  ice  salinity.  This  is  done  by  mul¬ 
tiplying  the  masses  of  all  the  components,  except 
pure  ice,  by  Sj/34.325  where  5(  is  the  ice  salinity  in 


parts  per  thousand.  The  mass  of  pure  ice  is  then  1000 
g  minus  the  total  ma>s  of  brine  and  precipitated  salts. 
Even  though  these  calculations  are  straightforward, 
they  are  both  time-consuming  and  prone  to  arithmetic 
errors.  In  fact  they  are  rarely  performed  and  estimates 
of  the  amount  of  gas  present  in  sea  ice  and  of  the 
total  void  volume  of  the  ice  are  rarely  obtained. 

In  the  present  report,  equations  are  derived  from 
which  the  gas  volume  and  brine  volume  can  be  di¬ 
rectly  calculated  given  the  ice  salinity,  temperature, 
and  density.  The  equations  cover  the  temperature 
range  from  -2°C  to  -30°C  and  consider  the  presence 
of  solid  salts.  Equations  are  also  given  to  calculate 
sea  ice  density  and  the  variation  of  gas  and  brine  vol¬ 
ume  of  a  sea  ice  sample,  if  the  temperature  of  the 
sample  is  changed  from  the  temperature  at  which  the 
ice  density  was  originally  obtained. 

PREVIOUS  WORK 

Only  a  few  investigators  have  advanced  methods 
for  calculating  the  density  and  gas  volume  in  sea  ice. 
Zubov  (1945)  calculated  the  density  of  gas-free  sea 
ice  and  provided  a  table  of  values  that  could  be  used 
at  temperatures  of  -23°C  and  higher.  However,  his 
work  was  done  prior  to  Nelson  and  Thompson’s 
(1954)  experiments  on  the  chemistry  of  sea  ice  brine. 
Zubov  also  did  not  take  into  account  the  presence  of 
solid  salts.  Anderson  (1960)  published  a  table  giving 
the  density  of  gas-free  sea  ice  at  14  selected  tempera¬ 
tures  in  the  temperature  range  between  0  and  -50°C. 
This  work  was  based  on  Assur’s  (1958)  phase  relations. 
As  the  exact  method  of  calculation  is  not  given,  it  is 
not  clear  whether  Anderson  considered  the  presence  of 


solid  salts  in  the  ice.  The  only  relation  available 
that  calculates  the  gas  volume  of  sea  ice  directly  was 
developed  by  Schwerdtfeger  ( 1 963).  However,  it  is 
only  valid  at  temperatures  warmer  than  -8.2°C,  the 
Na2S04  •  10  HjO  crystallization  temperature.  He 
also  neglected  the  presence  of  solid  salts  and  assumed 
that  the  volume  of  brine  was  equal  to  that  of  pure 
water.  For  warm  sea  ice  these  are  reasonable  approx¬ 
imations.  At  the  conclusion  of  this  report  compari¬ 
sons  are  made  between  the  results  of  these  previous 
authors  and  the  present  results. 

It  should  also  be  noted  here  that  Assur's  (1958) 
brine  volume  table  and  the  brine  volume  equations 
based  on  this  table  derived  by  Frankenstein  and 
Garner  (1967)  assume  a  constant  ice  density  of  0.926 
Mg/m3.  Consequently,  these  reports  should  be  used 
with  caution  in  gas  volume  calculations.  The  brine 
volume  obtained  from  these  sources  should  be  mul¬ 
tiplied  by  p/0.926,  where  p  is  the  gas-free  theoretical 
density,  or  better  still,  the  measured  ice  density  in 
Mg/m3  at  the  temperature  of  interest,  to  determine 
the  brine  volume  for  a  given  sea  ice  density. 


DERIVATION  OF  EQUATIONS 

A  sample  of  sea  ice  is  composed  of  pure  ice,  brine, 
solid  salts,  and  gas.  In  determining  the  gas  volume 
of  sea  ice  we  are  interested  in  the  relative  masses  and 
volumes  of  the  different  components  at  a  given  ice 
salinity,  temperature,  and  bulk  density. 

In  the  derivations  to  follow,  mt,  p8  and  Vt  are 
the  mass,  density,  and  volume  of  component  8  where 
the  subscripts  a,  b,  /',  s,  and  ss  denote  the  component 
air,  brine,  pure  ice,  salt,  and  solid  salts,  respectively. 
M  is  the  bulk  mass,  V  the  bulk  volume,  and  p  the 
bulk  density.  These  should  not  be  confused  with  the 
gas-free  theoretical  mass,  volume,  and  density.  The 
terms  mb  and  m5*  denote  the  mass  of  salt  in  the  brine 
and  mass  of  salt  in  the  solid  salts,  respectively.  The 
mass  of  gas  in  the  ice  is  assumed  to  be  negligible. 

Brine  volume 

The  salinity  of  the  iceSj  is  defined  as 
mb  +  m\s 

c  -  _L=  - : - 1 - 

'  M  mb  +ms5  +Z77j 

The  salinity  of  the  brine  Sb  is  defined  as 
mb 

From  these  two  relations  an  equation  for  the 
brine  volume  of  sea  ice  can  be  derived.  It  should  be 


(D 

(2) 


noted  that  the  brine  salinity  and  the  relative  amounts 
of  salt  in  the  brine  and  solid  salts  are  unique  functions 
of  ice  temperature  via  phase  relations. 

From  eq  1 


MS  j  =  +  mb 


and 


or 


mb  =  MS  j  -  m5* 


mb  =  MS |  -  kmb 


where 


(3) 


Solving  for  mb  from  eq  3  and  noting  from  eq  2 
that 

"K  =  Pb 

we  obtain 


PS  i 

^  =  Pb^b 


fe) 


which  because  M  =  pV can  also  be  written  as 


pS, 


V  Pb5b 


(4) 


where  VbfV  is  the  relative  brine  volume.  If  we  neglect 
the  salt  present  as  solid  salt,  k  =  0,  eq  4  reduces  to 


^  _?£l_ 

*  '  Pb*b  ’ 

a  relation  utilized  earlier  by  Weeks  and  Lofgren  (1967). 

Because  the  brine  density  can  be  approximated 
(Zubov  1945,  and  Cox  and  Weeks  1975  by 

pb  (Mg/m3)  =  1  +  0.0008  Sb  (<fao) 

and  because  5b  and  k  are  unique  functions  of  temper¬ 
ature,  the  relative  brine  volume  can  be  expressed  as 


v*  f£i_ 

V  “  F,(D 

where  5j  is  in  °/oo,  p  is  in  Mg/m3,  and 


(5) 


Table  1.  Values  foe  Sb,  k,  C,  F,  { T ),  and  F2(7)  at 
different  temperatures. 


F\(T)~  Pb5b  (1+*). 


(6) 


T 

rej 

s* 

<°foo) 

k  X  103 

CX  103 

(M|/m3) 

f2Ct) 

-  2 

37.6 

0 

0 

38.731 

0.123 

-  4 

70.6 

0354 

0.148 

74.662 

0.151 

-  6 

99.8 

1.050 

0.387 

107.876 

0.177 

-  8 

126.5 

1.400 

0.660 

139.441 

0.199 

-to 

142.8 

55.277 

18.256 

167.865 

0.222 

-12 

IS  7.6 

84.141 

30.493 

192.378 

0.240 

-14 

1713 

97.627 

38.421 

214.143 

0.2S6 

-16 

184.4 

106.330 

44.952 

234.033 

0.271 

-18 

197.0 

112370 

50.808 

2S  3.588 

0.284 

-20 

209.9 

118.078 

56.8S1 

274.074 

0.298 

-22 

222.6 

123.090 

63.015 

294.496 

0.312 

-24 

230.5 

509.787 

217.168 

412.236 

0.394 

-26 

232.7 

1312.694 

537.697 

638.433 

0.556 

-28 

234.1 

2065.827 

842.341 

852.171 

0.708 

-30 

235.6 

2685.708 

1098.887 

1032.102 

0.836 

Values  for  F^(T)  are  given  in  Table  I  and  plotted 
against  temperature  in  Figure  1.  Sb  and  k  are  de¬ 
termined  from  Assur’s  (1958)  phase  equilibrium 
table.  The  lines  throu^i  the  data  in  Figure  1  are 
least-squares  curves  of  the  form 

F\(T)  =  a0  +  ot,7'  +  a2  T2  +  a37‘3  (7) 

where  T  is  the  ice  temperature  (°C).  Coefficients 
for  the  curves  are  given  in  Table  2. 

Solid  salt  volume 

Also,  via  the  phase  relations,  at  any  given  temper¬ 
ature  the  mass  of  solid  salts  (if  solid  salts  are  present) 
is  proportional  to  the  mass  of  brine;  that  is 

=  Cmb  (g) 


Figure  1.  Function  F  y(T)  versus  temperature.  Curve  determined  by  method 
of  least-squares. 


Table  2.  Coefficients  for  functions  F,  (T)  and  F2(7}  determined  by  method  of 
least-squares. 


Correlation 


T  CC) 

“0 

a2 

a3 

coefficient 

Fy  (71 

-2>7>-22,9 

-4.732 

-2.245  X 

101  -6.397  X 

10"1 

-1.074  X  10*2 

0.9999 

-22.9>7>-30 

9.899  X  103 

1.309  X 

103  5.527  X 

101 

7.160  X  10*1 

0.9999 

F2(71 

-2>7>-22.9 

8.903  X  10-2 

-1.763  X 

10-2  -5.330  X 

10-4 

-8.801  X  10*6 

0.9999 

-22.9>7>-30 

8.547 

1.089 

4.518  X 

10‘2 

5.819  X  10-* 

0.9999 
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where  C  is  a  function  of  T,  The  relative  volume  of 
solid  salts,  VK/  V,  is  then  equal  to 


V 


(9) 


where  p$s  is  the  average  solid  salt  density  assumed  to 
be  constant  at  1.5  Mg/m3. 


Pure  ice  volume 

The  mass  of  pure  ice  is  equal  to 


m,  =  M  -  mb  -  mis 
or  from  eq  8 

m.=M  -  (1+C)/wb  . 


Since 


m\  =  P\  V\ 


M  =  pV 


mb  =PbVb 


we  have 


p^pV-O+C)  PbVb. 


Solving  for  the  relative  pure  ice  volume,  VJV,  we 
obtain: 


Zlm£-.h+ci  ^  ^ 

V'  Pi  '  '  Pi  v 


(10) 


Air  volume 

The  relative  air  volume,  VJV,  is  equal  to 


Substituting  eq  9  and  1 0  into  1 1 ,  we  obtain 


(12) 


Equation  12  can  be  simplified  by  defining 

'w|',4°  ?  l,3) 

substituting  for  (VJV)  from  eq  4,  and  recalling  the 
definition  of  F,  (T)  given  in  eq  6.  The  resulting  re¬ 
lation  is 


li  =1  .£- 
V  p , 


pSj  F2(T) 
FJT) 


(14) 


Values  for  F2(T)  are  given  in  Table  1  and  plotted 
against  temperature  in  Figure  2.  C  is  determined 
from  Assur’s  (1958)  phase  equilibrium  table.  The 
lines  through  the  data  in  Figure  2  are  also  least-squares 
of  the  same  form  as  eq  7. 

Coefficients  for  these  curves  are  given  in  Table  2. 

In  determining  the  value  of  F2(T)  at  a  given  tempera¬ 
ture,  the  pure  ice  density  was  calculated  from  (Pounder 
1965): 

Pi  (Mg/m3)  =  0.917  -  1.403  x  10-4F(°C) . 


Figure  Z  Function  F2(T)  vs  temperature.  Curve  determined  by  method  of 
least-squares. 


Sea  ice  density 

The  sea  ice  density  can  be  found  from  eq  14  by 
solving  for  p: 

L  Pi  W  .... 

p~  1  v)  •  (15) 

Unlike  the  calculated  densities  given  in  Schwerdtfeger 
(1963)  and  Anderson  (I960),  this  equation  consid¬ 
ers  the  presence  of  both  air  and  solid  salts  in  the  ice. 

A  comparison  between  the  different  density  estimates 
is  given  later  in  this  report. 


From  eq  18  we  can  also  conclude  that 


Vb  lV-Vt)F3Qu7) 

V^=  (V*  ~  k'j)  F\  (S,,  7)  • 


(20) 


Substituting  eq  19  into  eq  20  and  solving  for  VJV 
we  obtain  for  the  new  relative  air  volume 


/.  K\,pJ^(S,T)F\  (T) 

V  V  Vrp')\F3(Si,T)F,  (T) 

(21) 


Gas  volume  at  another  temperature 

The  gas  volume  and  brine  volume  of  sea  ice  are 
determined  by  measuring  the  volume  and  mass  of 
the  ice  sample  at  a  given  temperature  and  the  salin¬ 
ity  of  the  melted  sample.  If  tests  are  conducted  at 
a  temperature  different  than  the  temperature  at 
which  the  bulk  density  was  determined,  the  air  and 
brine  volumes  must  be  recalculated. 

Equations  are  now  derived  to  calculate  the  air 
and  brine  volumes  at  a  temperature  that  differs  from 
the  temperature  at  which  the  ice  density  was  de¬ 
termined.  Primed  variables  denote  the  initial  value 
of  the  different  variables. 

It  is  initially  assumed  that  the  brine  and  air  pock¬ 
ets  are  interconnected.  From  eq  5  we  know  that 

A#  S= 

Vb=F^n  ■ 


Because  the  mass  and  bulk  salinity  of  the  ice  do  not 
change  with  temperature 


F\  (T) 


V*  ^  (T)  ■ 


(16) 


By  rearranging  eq  14,  substituting  for  p  based  on  eq 
5  and  defining 


F3(sit-n  = 


pi5i 

Fycn-p^F^n 


(17) 


it  can  be  shown  that 


Vb*(V-Vt)F3$x,r\.  (18) 


Also  dividing  both  sides  of  eq  16  by  V  and  noting  that 


P'V 

P 


If  the  volume  of  the  sample  is  not  known  at  the  new 
temperature,  one  can  assume  that  the  change  in  the 
external  dimensions  of  the  sample  is  due  to  the  vol¬ 
umetric  change  in  the  ice  matrix;  that  is 


In  the  temperature  range  where  many  investigations 
are  carried  out 


Pi 

I 

Pi 


1  . 


If  the  brine  and  gas  pockets  are  not  connected,  a 
change  in  gas  volume  may  also  take  place  if  the  ice  is 
warmed,  that  is,  as  less  dense  ice  is  melted  to  dilute 
the  brine.  In  this  case  the  new  relative  air  volume  is 
equal  to 


^  K  ,  p  n(S.„T)F\  (T) 
V~  v+'~7  f3  (S,,  T)  F,  (7) 


(22) 


and  in  many  applications 


V-  -4*1. 

P  P; 

If  the  brine  and  air  pockets  are  not  connected  and 
the  ice  is  cooled,  the  relative  air  volume  will  not 
change.  However,  due  to  phase  changes,  some  brine 
will  be  lost  from  the  sample  as  a  result  of  brine  ex¬ 
pulsion.  The  new  relative  brine  volume  can  be  cal¬ 
culated  from  eq  19. 


DISCUSSION 


we  obtain  for  the  new  relative  brine  volume 


^  p  ^  F\  m 

V  ~  7  V  Fj( 7)  ■ 


(n) 


As  methods  for  calculating  the  gas  content  of  sea 
ice  invariably  contain  a  calculation  of  the  density  of 
gas-free  sea  ice,  it  is  interesting  to  compare  the  den¬ 
sity  estimates  developed  in  this  report  with  previously 
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Table  3.  The  calculated  densities  (Mg/m3)  ,f  gas-free  sea  ice 
at  different  temperatures  and  salinities  as  given  by  different 
authors. 


1 

3 

5 

10 

20°  loo 

1 

u 

V. 

1 

II 

Zubov  (1945) 

0.922 

0.926 

0.930 

0.939 

Anderson  (1960) 

0.9198 

0.9252 

0.9307 

0.9447 

0.9739 

Schwerdtfeger  (1963) 

0.9191 

0.9233 

0,9275 

0.9379 

0.9588 

This  report 

0.9200 

0.9254 

0.9308 

0.9448 

0.9740 

T  =  -8°C 

Zubov (1945) 

0.920 

0.923 

0.925 

0.932 

0.944 

Anderson  (1960) 

0.9192 

0.9216 

0.9240 

0.9301 

0.9426 

Schwerdtfeger  (1963) 

0.9175 

0.9186 

0.9196 

0.9222 

0.9275 

This  report 

0.9193 

0.9218 

0.9242 

0.9303 

0.9429 

T=  -10°C 

Zubov (1945) 

0.920 

0.922 

0.925 

0.929 

Anderson  (1960) 

0.9194 

0.9216 

0.9239 

0.9296 

0.94 1 2 

This  report 

0.9195 

0.9218 

0.9240 

0.9297 

0.9413 

r  =  -30°c 

Anderson  (1960) 

0.9216 

0.9230 

0.9245 

0.9281 

0.9354 

This  report 

0.9219 

0.9233 

0."247 

0.9281 

0.9352 

published  estimates.  Table  3  gives  such  a  compari¬ 
son.  As  can  be  seen,  the  agreement  between  the 
present  results  and  those  of  Anderson  (1960)  is  very 
good  (the  present  values  are  at  most  +0.0003  Mg/m3 
higher  with  a  modal  difference  of  only  +0.0001 
Mg/m3).  This  agreement  is,  however,  not  surprising 
in  that  both  Anderson’s  and  our  calculations  are 
based  on  Assur’s  phase  equilibrium  tables.  The 
agreement  essentially  verifies  both  sets  of  calcula¬ 
tions.  When  comparisons  are  made  with  Schwerdt- 
feger’s  (1963)  results  for  sea  ice  warmer  than  -8.2°C, 
much  larger  differences  are  found.  Within  this  tem¬ 
perature  range  these  differences  vary  from  +0.0009 
at  warm  temperatures  (-2°C)  and  low  salinities 
(1  °/oo)  to  +0.0154  Mg/m3  at  -8°C  and  20  °/oo.  The 
agreement  with  the  earlier  values  of  Zubov  (1945) 
is  better  than  with  Schwerdtfeger,  with  the  largest 
difference  being  +0.0058  Mg/m3. 

It  is  also  interesting  to  use  the  relations  that  have 
been  developed  to  examine  the  importance  of  con¬ 
sidering  the  presence  of  solid  salts  in  making  gas  vol¬ 
ume  calculations.  For  sea  ice  warmer  than  -8.2°C, 
Schwerdtfeger  (1963)  calculated  the  relative  air  vol¬ 
ume  from 


V.  o  /1000-Sj  4.98  Si  \ 

~V  ~  1  ‘  1000  (  6.917  +  T  /  ( 2 3* 

where  p  is  in  Mg/m3,  S|  in  °/oo,  and  T  in  °C.  As 
mentioned  earlier  he  neglected  the  presence  of  solid 
salts  and  assumes  that  the  volume  of  brine  in  the  ice 
is  equal  to  the  volume  of  pure  water.  If  we  neglect 
the  presence  of  solid  salts  in  our  equations,  eq  12 
simplifies  to 


The  results  from  eq  23  and  24  arc  compared  to  the 
results  from  eq  1 2  in  Table  4.  The  air  volume  of  sea 
ice  of  different  assumed  salinities,  densities,  and  tem¬ 
peratures  are  calculated  using  these  three  equations. 
Schwerdtfeger's  results  are  given  in  column  1 ,  our 
results  neglecting  solid  salts  in  column  2,  and  our 
results  considering  solid  salts  in  column  3. 

All  the  calculated  values  at  a  given  salinity,  density, 
and  temperature  are  in  close  agreement.  This  indi¬ 
cates  that,  even  for  sea  ice  below  the  NaCI-2H20 
eutectic  temperature,  neglecting  the  presence  of  solid 


(' 
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Table  4.  Air  volume  (^oo)  calculations  from  three 
methods  for  sea  ice  having  difference  densities,  salin¬ 
ities  and  temperatures  (see  text  for  explanation). 

_ _ pfMg/m3) _ 


T  °'890  a9w  0-930 


Cc) 

1 

2 

3  1  2 

3 

1 

7 

7 

Sj  =  1  °A» 

-  6 

31.1 

31.8 

31.8  9.4  10.0 

10.0 

_ 

. 

-10 

- 

32.0 

32.1  -  10.2 

10.4 

- 

_ 

_ 

-20 

- 

33.4 

33.4  -  11.7 

11.7 

_ 

_ 

-30 

- 

34.8 

34.6  -  13.1 

12.9 

- 

- 

- 

S;  =  10  %o 

-  6 

46.5 

44.9 

44.8  25.1  23.5 

23.3 

3.7 

2.0 

1.9 

-10 

- 

42.7 

42.7  -  21.2 

21.2 

_ 

_ 

-20 

- 

42.2 

42.1  -  20.7 

20.6 

_ 

_ 

-30 

- 

43.1 

41.0  -  21.6 

19.4 

- 

- 

salts  in  air  volume  calculations  is  a  reasonable  assump¬ 
tion.  For  warm,  high  density,  high  salinity  sea  ice 
(0.930  Mg/m3),  Schwerdtfeger’s  value  differs  by  1.8 
^oo,  or  nearly  100%.  However,  the  brine  volume  of 
this  ice  is  relatively  much  greater  (85.1  ^oo)  and 
there  is  little  difference  in  the  total  ice  porosity. 

CONCLUSIONS 

Equations  have  been  derived  to  quickly  calculate 
the  air  volume  of  sea  ice  considering  the  presence 
and  absence  of  solid  salts.  The  results  show  that 
neglecting  the  presence  of  solid  salts  is  a  reasonable 
assumption.  However,  in  this  era  of  computers  and 
sophisticated  hand  calculators,  we  recommend  con¬ 
sidering  the  presence  of  solid  salts  and  the  use  of  eq 
12  to  calculate  the  air  volume  of  sea  ice.  A  program 
for  making  such  calculations  on  a  Hewlett-Packard 
41 C  calculator  is  given  in  Appendix  A. 

Earlier  we  mentioned  that  very  few  direct  deter¬ 
minations  of  the  gas  content  of  sea  ice  have  been 
undertaken.  This  might  lead  the  reader  to  question 
the  usefulness  of  the  calculations  outlined  here  be¬ 
cause  of  a  lack  of  available  comparisons  between  cal¬ 
culated  gas  porosities  and  measured  gas  porosities. 
Fortunately  such  a  comparison  has  now  been  com¬ 
pleted  for  first-year  sea  ice  (Nakawo  1983)  indicating 
a  very  good  agreement  between  the  calculated  and 
measured  gas  porosities.  Inasmuch  as  the  calculation 
procedure  we  have  developed  is  very  rapid  and  the 
measurements  required  to  utilize  it  are  simple,  it  is 


hoped  that  gas  volume  will,  in  the  future,  join  brine 
volume  and  specimen  temperature  as  part  of  the 
sample  description  for  all  physical  property  tests  on 
sea  ice. 
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APPENDIX  A.  HP-41  C  AIR  VOLUME  PROGRAM 


The  following  HP-41C  program  was  written  to  calculate  the  gas  volume  ( VA),  brine  volume  (VB), 
and  total  porosity  ( N )  of  sea  ice  samples.  The  program  firs:  prompts  the  user  for  the  ice  density 
(g/cm3),  salinity  (°/oo),  and  temperature  (°C)  at  which  me  ice  density  was  determined.  The  gas  vol¬ 
ume,  brine  volume,  and  total  porosity  are  then  calculated  for  this  temperature.  Next  the  program 
prompts  the  user  for  the  test  temperature  and  calculates  the  test  gas  volume,  brine  volume,  total 
porosity,  and  new  density.  It  is  assumed  that  the  brine  and  air  pockets  are  interconnected. 

To  initialize  the  program,  the  contents  of  storage  registers  32  to  50  are  entered.  The  program 
is  then  executed  and  after  each  real  time  input  (PROMPT),  RUN  is  entered.  During  execution, 
storage  registers  1  through  15  are  used.  A  sample  jtput  is  also  given  below. 


MAIN  PROGRAM 


PRP  -flIRVOL* 

38  • 

80»LBL  03 

39  RCL  04 

81  XEQ  13 

81*L8L  -filPVCL* 

40  / 

82  STO  09 

02  *1IEN  G/CH- 

41  ♦ 

83  XEQ  14 

03  PROMPT 

42  STO  07 

84  STO  18 

04  STO  81 

43  RCL  58 

05  -SOL  PPT  ?■ 

44  * 

85*LBL  84 

06  PROMPT 

45  FIX  1 

86  RCL  06 

07  STO  02 

46  *VQ  PPT  =  * 

87  RCL  02 

08  'TEMP  C  ?■ 

47  RRCl  X 

80  • 

09  PROMPT 

48  RVIEN 

89  RCL  04 

10  STO  03 

49  RCL  01 

90  ENTERt 

11  -22.9 

50  RCL  82 

91  RCL  06 

12  XOY 

51  • 

92  RCL  02 

13  XOY? 

52  RCL  04 

93  » 

14  GTO  01 

53  / 

94  RCL  85 

15  XEQ  11 

54  STO  88 

95  * 

16  STO  04 

55  RCL  50 

96  - 

17  XEQ  12 

56  • 

97  / 

18  STO  05 

57  *VB  PPT  =  * 

98  STO  11 

19  GTO  82 

58  RRCL  X 

99  XEQ  15 

59  RVIEN 

100  STO  12 

28*LBL  81 

60  RCL  07 

101  RCL  02 

21  XEQ  13 

61  RCL  08 

182  * 

22  STO  04 

62  ♦ 

103  RCL  09 

23  XEQ  14 

63  RCL  50 

104  ENTERt 

24  STO  05 

64  * 

105  RCL  12 

65  *H  PPT  *  • 

106  RCL  02 

25«LBL  02 

66  RRCL  X 

107  * 

26  XEQ  15 

67  RVIEN 

108  RCL  10 

27  STO  06 

68  'TEST  TEMP  C  ?* 

189  * 

28  1 

69  PROMPT 

110  - 

29  ENTERt 

78  STO  03 

111  / 

30  RCL  01 

71  -22.9 

112  STO  13 

31  RCL  06 

72  XOY 

113  1 

32  / 

73  XOY? 

114  ENTERt 

33  - 

74  GTO  03 

115  1 

34  RCL  01 

75  XEQ  11 

116  ENTERt 

35  RCL  02 

76  STO  09 

117  RCL  07 

36  * 

77  XEQ  12 

118  - 

37  RCL  05 

78  STO  18 

119  RCL  11 

79  GTO  04 

120  * 
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macixatta  pads  bummwt  filmed 


121  RCL  04 

122  * 

123  RCL  13 

124  / 

125  RCL  09 

126  / 

12?  RCL  12 
128  ♦ 

129  RCL  06 
138  / 

131  - 

132  STO  14 

133  RCL  58 

134  * 

135  'TEST  Vfl  PPT  = 

136  ARCL  Y 

137  AVIEM 

138  1 

139  ENTEPt 
148  RCL  14 

141  - 

142  RCL  13 

143  * 

144  STO  15 

145  RCL  58 

146  * 

147  'TEST  VB  PPT  =  ■ 

148  ARCL  X 

149  AVIEH 

150  RCL  14 

151  RCL  15 

152  + 

153  RCL  58 

154  * 

155  'TEST  N  PPT  =  ■ 

156  ARCL  X 

157  AVIEK 

158  RCL  15 

159  RCL  89 

160  « 

161  RCL  02 

162  / 

163  FIX  3 

164  ’TEST  DEN  C'CC  =' 
165  ARCL  X 


166  flVIEN 
167  CTO  'AIRV9L- 
168«LBL  11 
169  RCL  32 
178  ENTERT 

171  RCL  03 

172  RCL  33 

173  * 

174  ♦ 

175  RCL  03 

176  Xt2 

177  RCL  34 

178  * 

179  ♦ 

180  RCL  03 

181  3 

182  YTX 

183  RCL  35 

184  • 

185  ♦ 

186  RTN 

187*LBL  12 

188  RCL  36 

189  ENTERt 
196  RCL  37 

191  RCL  83 

192  « 

193  ♦ 

194  RCL  03 

195  Xt2 

196  RCL  38 

197  * 

198  ♦ 

199  RCL  03 

280  3 

281  YTX 

202  RCL  39 

203  • 

204  ♦ 

205  RTN 

286*LBL  13 

207  RCL  40 

208  ENTERT 


289  RCL  41 
218  RCL  13 
211  • 

212  ♦ 

213  RCL  03 

214  XT2 

215  RCL  42 

216  # 

217  ♦ 

218  RCL  83 

219  3 
228  YTX 

221  RCL  43 

222  • 

223  ♦ 

224  RTN 

225*L0L  1* 

226  RCL  44 

227  ENTERT 

228  RCL  45 

229  RCL  83 
238  « 

231  ♦ 

232  RCL  83 

233  XT2 

234  RCL  46 

235  • 

236  ♦ 

237  RCL  83 

238  3 

239  YTX 

240  RCL  47 

241  • 

242  ♦ 

243  RTN 

244*LBL  15 

245  RCL  48 

246  ENTERT 

247  RCL  49 

248  RCL  03 

249  * 

250  ♦ 

251  RTN 

252  END 


INITIALIZATION  OF  THE  PROGRAM  SAMPLE  RUN 


832.850  ENTEP‘  XEQ  *MW0L’ 


PRREGX 

BEN  G/CC? 

.919 

RUN 

R32= 

-4.732+00 

SAL  PPT  ? 

R33= 

-2.245+81 

4.588 

RUN 

R34= 

-6.397-01 

TEW  C  ? 

R35= 

-1.074-02 

-15.000 

RUN 

R36= 

8.983-02 

VP  PPT  =  14.7 

R37= 

-1.763-82 

VB  PPT  *  18.3 

R39= 

-5.338-04 

N  PPT  =  33.0 

R39= 

-8.801-86 

TEST  TEW  C  ? 

R40= 

9.899+03 

-5.8 

RUN 

R41= 

1.309+03 

TEST  VP  PPT  =  17.2 

R42= 

5.527+61 

TEST  VB  PPT  =  44.0 

R43= 

7.160-01 

TEST  N  PPT  =  61.2 

R44= 

8.547+00 

TEST  1EN  G/CC  =0.989 

R45= 

1.039+00 

BEN  G/CC? 

R46=  4.518-02 
R47=  5.819-04 
R4S=  9.170-01 
R49=  -1.403-04 
R58=  1.080+03 
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A  facsimile  catalog  card  in  Library  of  Congress  MARC 
format  is  reproduced  below. 
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